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A mathematical model of the waterjet propulsion 
plant for surface effect ships consisting of gas turbine, 
reduction gear, axial pump, variable flushed inlet, diffuser 
and fixed nozzle was developed. Control algorithms were 
also developed for the variable geometry flushed inlet and 
for the pump RPM, and were tested by digital simulation of 
the mathematical model. Pump cavitation was included in the 
model and simulation. 


The control scheme for the variable flushed inlet 
involves the maintainance of a diffusion flow pattern about 
the flushed inlet. The control scheme for the pump involves 
the maintainance of pump RPM as dictated by the throttle RPM 
setting. In both control schemes, the rate and hard limits, 
on the hydraulic actuator for the flushed inlet and on the 
fuel flow rate of the gas turbine, were considered. Results 
of real time tests of the developed algorithm are presented. 
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INTRODUCTION 


The inherent geometrical flexibility of Waterjet 
Propulsion system commands a significant interest for its 
use in Surface Effect Ship (S.E.S.). In particular where 
the power level of the propulsion train is high, waterjet 
system, despite its relatively lower propulsive efficiency 
compared to the controllable reversible pitch supercavi- 
tating propellers (Fig. 1) is considered to be a more 
reliable system to fit the sidewall configuration of the 
S.E.S. A great risk is involved when propeller shafts at 
power levels of 20,000 H.P. and above are turned around 
corners to drive propellers. 

Thus the current U.S.A. effort for the projected 2,000 
ton S.E.S. having a cruise speed of approximately 80 knots 
focuses attention to the waterjet system as the best choice 
for propulsion. As a result of the experience of the 100 
ton waterjet driven S.E.S., it was recognized that a vari- 
able inlet geometry was needed to maintain high pressure 
recovery at all speed levels, particularly at hump speed and 
at cruise (Fig. 2). Fig. 3 shows the typical speed resis- 


tance of the S.E.S. One of the design requirements is that 
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at hump speed the waterjet system must be able to deliver 

a maximum thrust of at least 20% greater than hump total 
resistance for an acceptable acceleration performance. 

This implies that the water flow rate for the system at 
cruise speed and hump speed will be approximately the same. 
Since the ratio of cruise speed to hump speed is approxi- 
mately 3, the flush inlet area must have to be adjusted 
accordingly so that a diffusion flow pattern will be main- 
tained always from the region in front of the flush inlet 
to the pump inlet (Fig. 4). If not, then pressure recovery 
will be seriously degraded and the pump may cavitate due to 
insufficient suction head pressure. 

At present, work is being conducted at M.I.T. to 
optimize the design of components of the waterjet system for 
the S.E.S. so maximum payload can be obtained for a given 
displacement range, cruise speed and hump speed. 

The purpose of this paper then is to model and simu- 
late the dynamic performance of these conponents and to 
identify and bring out problems that may arise in the overall 
control of the system and possible criteria for the selection 
of the components. The case study presented herein is that 
of a 2109 ton S.E.S. whose waterjet system consists of six 


LM2500 gas turbines, reduction gears, axial water pumps, 
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variable inlet geometry flushed inlet and fixed geometry 
nozzle. The simulation program has been set up in such a 
way that thrust reverser characteristics can also be added 


to analyze the backing performance of the craft. 





CHAPTER I 


WATERJET PROPULSION COMPONENT MA'‘TCHING 
I-1 GENERAL: 


In order for the waterjet system to function effi- 
ciently, pump design point must match accordingly to the 
thrust requirements both at hump and at cruise speeds. This 
includes the determination of the following: design range 
of variable inlet, area design pump RPM, design pump flow 
rate, design nozzle area. Using the momentum equation the 


thrust delivered by the system is: 


Tg = ee (V5 cosa - Vey) (1-1) 


where: Ty = thrust develop (1lbE£) 

fP = water density (lbm/£t3) 
= water flow rate (ft/sec) 
g = gravity. (32.2, £t/sec) 
V. = waterjet velocity (ft/sec) 
Vo = S.E.S. velocity (ft/sec) 
o& = depression angle of the nozzle 
For most design cosm is approximately equal to unity 


thus can be neglected in the preliminary calculation. The 
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thrust develop then is: 
QV, | V5 
ry = late] i 
The thrust required at hump speed is: 


T a (Roy + Ray) (1-3a) 


x. 
0) 
4K 
oO 
rH 
Il 


thrust required at hump 


RH speed (1b£) 

eer bare hull resistance at hump 
H speed (lbf) 

Rig = inlet appendage drag at hump 


speed (lbf) 


The factor 1.2 is the required minimum acceleration 
margin. 


fhe thrust required at cruise speed is: 


T = iP + R (1-3b) 


= 
=p 
@ 
R 
0) 
rH 
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thrust required at cruise 


Bb speed (lbf) 
BOAT bare hull resistance at cruise 
speed (lbf) 
Ryc = inlet appendage drag at cruise 


speed (lbf) 
The inlet appendage drag for equations (1-3) is: 


= to 
IH pH 2) YoHH (1-4a) 
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inlet drag coefficients 
at hump and cruise 
speeds respectively 


V and V = S.E.S. velocity at hump 
and cruise speeds 
respectively 


and Qa = water flow rate at hump 
and cruise speeds 
respectively 


I 


Equating thrust delivered and thrust required and 


substituting equations (1-4) to equations (1-3) results in: 


1.2 Rey 
QO, = 5 (1-5a) 
i VjH 
#0 Rsc 
Cae: Vic | 
g Voc [ Ee = (1+0 .5C,.) (1-5b) 


The pump head required at hump and cruise speeds from 


Bernoulli equation are respectively: 


2 2 
V V- 2gH 
OH HH + -9"wh 
Apa 7 6g vit) Soa % al Pee 
2 y 
H eae = SSS aul (1-6b) 
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overall efficiency of inlet 
and diffuser 
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height of the nozzle above the 
waterline (ft) 
The effective horsepower required (EHP = P Hp2) OG 


hump and cruise speeds from equations (1-5) and (1-6) are: 


V5H i 2gHwr, 
Se sO LN Vor Vou? ~ “Loa (1-7a) 
H ~ (2) (550) Van 


eee (aes) sieve 


(1-7b) 
a - (1+0.5 cy) | 
Oc 

Note that equations (1-5), (1-6), and (1-7) are all 
expressed as a function of jet velocity ratio (V5/Vo) - In 
the case of a fixed inlet area the jet velocity ratio and 


flow rate at hump and cruise speeds can be related to the 


nozzle jet area (A,) by: 
Se = (1-8) 


Fig. 5 shows a typical propulsive efficiency 
(RV,/EHP) versus jet velocity ratio (V/V) for different 


-values of loss coefficient Kp where total head loss 
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Pe K,V;/23- Maximum efficiency is observed to occur 
for most practical designs between V,/Vy of Lee tO ue? sO. 
Likewise equation (1-7b) can be differentiated with 
respect to (V/V) to obtain minimum EHP. However, the 
result may not be compatible to the requirement of equations 
(1-5a), (1-6a) and (1-7a) because of the fixed nozzle area 
constraint unless the total hump resistance (Roy + Roy) is 
lower than some critical value. For the case where the 
total hump speed resistance is greater than this critical 
value the approach in finding the optimum cross sectional 
area of the nozzle and jet velocity ratios at hump and 
cruise speeds is as follows: 

a. Choose some multiple integral number of existing 
prime movers with a given shaft horsepower 
(SHP) rating. 

b. Multiply the maximum intermittent SHP available 
by the reduction gear efficiency and pump 
efficiency to obtain the maximum intermittent 
total EHP available. 

c. Substitute in equation (1-7a) and solve the 
quadratic equation for Se Vi.) The “two 


OH 


values of (V5 Vow) represent the maximum and 


minimum values of (V5/ Von) for which equation 


(1-7a) will be satisfied. 
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d. Solve for a using equations (1l-5a) and (1-8) 
using values of Vin’ You within the range 
obtained in step c. 

e. Solve for V5/Voc and check equation (1-7b) 
if the total sustained EHP is within acceptable 
limits. 

f. Iterate steps a to e until optimum number of 
engines, Aa propulsive efficiency is obtained. 

For the case where the total hump speed resistance 


is less than the critical value as dictated by equations 


(1-7a), (1-7b) and (1-8) the procedure becomes easier since 


eee d ivantzOn 16 te t 1 f 
- 8(V5C/Voc) need only to be se o zero to solve for 


optimum (V5~/Voc) and the rest follows. 
I-2 PRELIMINARY PUMP DESIGN: 


The procedure stated herein applies primarily to 
determine the rated pump RPM, pump diameter, annulus area 
and other related pump characteristics based on designed flow 
coefficient (%) and head coefficient (), designed flow rate 
(Q) and designed pump head (H)) - By ee ete 


Cenk aes 
x n 
P= y = NID (1-9) 
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where: oe = water axial velocity (ft/sec) 

U = blade tip velocity (ft/sec) 

Q = water flow rate (ft/sec) 

A, = pump annulus area (et) 

N = pump RPM 

D = pump tip diameter (ft) 

The axial pump annulus area (a) is: 
Tp2 

Rete = d/o) (1-11) 


Where: d = pump hub diameter (ft). 


The design flow rate from equations (1-5) and the 
design pump head from equation (1-6) can be substituted 
in equations (1-9), (1-10) and (1-11) to solve for A,, N, D 


for a given value of hub to tip diameter ratio. 
I-3 REDUCTION GEAR RATIO: 


The purpose of the reduction gear is to match the 
point of best efficiency of the pump to the point of best 
efficiency of the prime mover. Engine performance maps are 


usually expressed as plots of break horse power versus RPM 
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for various levels of specific fuel consumption (Fig. 6 and 
Fig. 7). Performance maps are also expressed as plots of 
torque versus RPM for various levels of fuel consumption 
(Fig. 8 and Fig. 9). In either case the point of best 
efficiency corresponding to power level required in equations 
(1-7) plus additional power for frictional losses and pump 
efficiency can easily be pinpointed. Thus the reduction 

gear ratio is determined as the ratio of the prime mover 

RPM at this point of best efficiency and the designed pump 


RPM (N). 
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CHAPTER II 


MODELING OF S.E.S. AND WATERJET PROPULSION PLANT 


II-1 GENERAL: 


In order to simulate the interaction of the com- 
ponents of the waterjet propulsion plant with the motion 
of the S.E.S., it is necessary to obtain the mathematical 
model of such an installation. The model would then allow 
the development of the appropriate control algorithm with- 
out actually employing the hardware involved. To obtain 
such a model a thermodynamic approach was chosen. Fig. 10 
shows the power flow block diagram of the waterjet pro- 
pulsion train. At this point it is recammended for the 
reader to take a look on ref. (1) and ref. (2) for the 
bond graft technique which will be used frequently in this 


paper. 


Cr 
TURBINE | Wy 





¢ |FLUSHED INLETL "“@p | WATERJET 
RESISTANCE | V, AND NOZZLE Q PUMP 


Figure 10. Block diagram of waterjet propulsion train 
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where: W = fuel flow rate demand (lbs/hr) 
W. = actual fuel flow rate (lbs/hr) 
Tr = power turbine torque (ftlb) 


W = power turbine speed (RPM) 


et, = pump torque (ftlbs) 
Ww, = pump speed (RPM) 
a = pump total head (ft) 


T = thrust (lbs) 


V. = S.E.S. velocity (ft/sec) 


The above diagram is interesting from the point of 
view of thermodynamic accounting of energies involved. The 
light broken arrows represent signal power level and the 
darkened arrows represent high power level transmission. 
The product of the variables affixed to the arrows multiplied 
by some constant represents power. For simplicity, the 
energy storage and energy dissipator elements associated 
with each component and power transmission are not included 
in the detailed modeling of each component which follows. 
Notice that all the energy associated with the fuel will 
in the end be dissipated into the environment primarily due 
to the S.E.S. resistance and secondarily but not avoidably 
due to the "resistances" inherent with the inefficiency of 


each component and inefficiency of power transmission. Thus 
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from the thermodynamic point of view, for the model to be 
valid in the steady state the power aut for each com- 
ponent must be greater than the power output in either 
direction. Likewise the first law and second law of 
thermodynamics must be satisfied during the transient 


conditions. 


The purpose of fuel is to provide energy primarily 
in the form of heat to the driving engine which is a gas 
turbine in this case study. The energy associated with the 
fuel can be expressed as heating value per unit weight (h). 


Fuel reaction with air mainly is as follows: 


T 
FP 
+ + + + 
2 CoHig 21 o5 Pa 12 ACO 18 H50 6 °, ont =a 
12 CO, + 18 HO + Q) + Q, (2-1) 
where: Tpp = flush point temperature 
Qy: Q, = heat evolve in the combustion 


process respectively 


The degree on which the above combustion process proceeds 
to the right is dependent on combustion temperature, fuel-air 
ratio and time on which the process is allowed to proceed in 


the combustion chamber of the gas turbine. 
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Thus the total heat evolved per unit fuel weight as 
function of time (h(t)) can be assumed in block diagram 


to be: 





where: Ga d= heating values per unit weight 
of fuel from equation (2-1) 


h(t) = total heating value of fuel per 
weight ' 
We = fuel flow rate 


C11 Ts = the exponential time constants of 
the above combustion process. 
cy and T. are very small terms but significant when con- 
vection time from fuel spray area of the combustion chamber 


to the turbine inlet is considered (Fig. 16). 
II-3 FUEL PUMP: 


The purpose of the fuel pump is to drive the fuel 
from the fuel tank into the combustion chamber of the gas 
turbine which is at a higher pressure. For the case of a 
positive displacement fuel pump (Fig. 9) represents the 


power flow diagram. 
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Figure 11. Pump power flow diagram 


where: (a input torque into the pump 


O = pump speed 
R, = pump mechanical frictional resistance 
roe = total lead developed by pump 
We = weight fuel flow rate 
Rp = fluid resistance associated with fuel 
value and fuel pipings 
C. = pump torque 
I, =fuel inertance associated in the fuel pipings 


m = modulo relating cy to Hop or modulo relating 


W@W to We 


TF = ideal transformer of mechanical power to 
fluid power 


1-junction = point of same © or the same We. 
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Since I, and Th are not independent for positive 
displacement pumps, I,, can be referred to the fluid side 


(or I, can be referred to the mechanical side). Similarly 


Ra can be referred to the fluid side. For the case where: 


a mW. = or H =m 2-2a 
te) aE op ie ( ) 
(b) Rn = bW® (2-2b) 
) . “ 
Ro = W + 2-2 
(c e° 2gC, (P Ao) ( c) 
then 
92-3) 
I aw 
2 fe iB 2 iL 
foo) et - - W 
(m i Pp) qt (m@ m bW, wel I (C)+ GC5 PRo)2)) 
where: b = mechanical frictional resistance 
coefficient 
S6= fluid frictional resistance coefficient 


in pipings 
A_ = orifice area of the fuel valve 


Co = orifice coefficient of the valve 


II-4 GAS TURBINE: 


Fig. 12 shows a schematic diagram of a typical air- 
craft gas turbine derivative. It is composed of two major 
units - the gas generator and the power turbine. The gas 


generator is a twin spool engine, consisting of low pressure 
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FIGURE [2 SCHEMATIC DIAGRAM OF FT-4 GAS TURBINE 
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compressor (L.P.C.), high pressure compressor (H.P.C.), 
combustion section, high pressure Ren (L.H.T.) and 
low pressure turbine. The power turbine (P.T.) is adapted 
to the gas generator by an annular diffusing transition 
duct. Thus the operating speed of the power turbine is 
independent of the gas generator speeds. The engine also 
comes with a fuel control package to ensure that starting 
and operating schedules of the gas turbines are met auto- 
matically without the gas generator going into surge and 
without subjecting the turbine inlet to excessive gas 
temperature. 

In order to coordinate the gas turbine package to the 
waterjet propulsion system, a mathematical model is developed 
as a part of this paper. 

For a satisfactory model, it can be assumed that the 
main accumulator of energy occurs in the gas pipe volumes 
and in the rotating parts. Although some heat energy is 
stored in the metal parts, its accumulation and discharge is 
carried out sufficiently slower. Thus during sharp load 
changes the transient is mainly from the dynamic interaction 
of rotors, volume and power turbine speed regulators. 

When a sharp load torque decrease occurs, the power 


turbine speed regulator initial action is to decrease quickly 
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the fuel flow to low value open up discharge value to 
release air intensively from gas See ae cer to evacuate 
the large quantity of stored energy. After which the 
discharge valve is closed and fuel flow is increased to 
the desired RPM. 

When an increase load is desired on the power turbine, 
it is necessary to increase the gas generator speed and 
therefore the stored energy therein. However in this 
transient process the compressor speed lags behind the 
pressure and temperature increase due to fuel flow increase. 
Thus not only for the protection of turbine against ex- 
cessive temperature but also preventing compressor insurge, 
the rate of fuel flow increase is strictly limited in 
accordance with air flow increase. Surge is also exciuded 
at low speeds by discharging air at the first group of stages 
or from compressor delivery valve. 

In order to see the basic principle underlying surge, 
consider a multistage axial compressor having the same flow 
coefficient @ = C/U and air axial velocity Ce for all 
stages, at design point. At very low speed ce has to increase 
greatly through the compressor since the outlet air density 
is now approximately equal to inlet air density. This means 


that the first stages are operating at a very high incidence 
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angle or positive stall value while the last stages will 

be at negative incidence angle or schist stall value or 
turbining value. Now, with sudden increase of fuel flow 
the compressor outlet pressure will also increase very 
quickly. Thus SS also decreases shifting the positive stall 
region toward the last stages. Since the first stages are 
not contributing any pressure increase, back air flow 
rapidly occurs as soon as stall hits the last stage with an 
accompanying explosive sound. 

With the above qualitative background, Fig. 13 is 
shown as the bond graph model of the gas turbine unit. 
Three variables are needed to describe power flow of com- 
pressible gas: mass flow rate (m.). total pressure (Piya 
total temperature (TO). The subscripts are evident from 


Fig. 12. The rest of the notations are: 


a 6 


I! 


torque 

© = angular speed 

We = fuel weight flow rate 

R = frictional resistance for shafts or orifices 
I = polar moment of inertia 

H_ = total head of fuel 


V = effective air volume capacitance 
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From the diagram the following are assumed: 


Pies as = ea on. es’ aa 
TGA fz TB! es ey aaa = See 27 on 


ola aol Loghe Hole Gato Pam 


Above assumption implies that pressure drops along ducts 


and diffuser sections of the gas turbines must be included 


in the overall efficiency of the adjacent compressors and 


turbines. 


From Newton's Law: 


dW;, 


Imy, dt = (CA = ont b,W,) 


dW 
Ing at = (Cap = Tue — }y Wy) 
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Assuming air follows ideal gas law: 


Pov = m RT, 


then: 





(2-4a) 


(2-4b) 


(2-4c) 


(2-5) 


(2-5a) 


(2-5b) 
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Vo itis (m,. - mg ™5) RT 4 Be a (2-5c) 
where: 
mia = CaPoa Y2g (Poa - Po.) (2-5d) 
dP 56 7 
3 at (ths - m6) RTs + mT. (2-5e) 
where: 
Mas = CgPAcgllos - Poc)29 
dm3 i, } 
dit a je ss 6 (2-5£) 
GP 67 . ‘ ° 
4 dt Es a a mg) QT O38 i We (2-5g) 
ane i 
dey Cay tas (2-Sh) 
dP 610 A ‘ A ; 
Me Pose Msg ~ Maj9) RTK19 + MERTK10 (2-5i) 
dm, 
ae" as Baw (2-53) 


Equations (2-4) and (2-5) are the equations governing 
the dynamics of the gas turbine in Fig. 12. 


From compressor and turbine maps Po Pol versus 


(31 J To1)/(ApP oz) for different reais co hea polytropic 
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efficiency are plotted (Fig. 14 and Fig. 15) Thus 
az, can be found as functions Poae 61" Tol and annulus 
area (A,). 

From Euler's equation, the torque on the rotor of 


an ideal lossless compressor or turbine is 


: mr ( &Ce) 
g 
In terms of power: 
= mo rAC, = The Aho 
Therefore: 
v a Mo1 Sho eae 
Lc > Wy, ( a 


Equation (2-6) will be true for an adiabatic compressor 
having no storage elements. Thus for a compressor having 
a polytropic efficiency 


( ¥ - 1) R 


—_— _ Cp 
™ polyc (n - 1) G5 
n n 
therefore: 
(i = 1). R 
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where 


To2 
Po = Tor (2-7) 


(2-8) 


Likewise for adiabatic turbines having no storage 
elements and with a give F 
i g n polyp 


(n-1) 
n 


— i a Oe ple de 
™ polyp -_ wists Ons n 





therefore 


Pp 
se) Cp 


-Ti = m Car Pag 


a8 p o8 


Plpoly » 
( | 


(2-9) 


Equations similar to (2-8) and (2-9) can likewise be 


developed for the H.P.C., H.P.T. and P.T., thus defining 


the differential equations (2-4). Since from equation (2-7) 


the total temperatures and from performance maps mass flow 








4¢ on vrekved Fook Ghd agseeiat hse 308 


; ray neeeey sp ite 
wrakeas 


LL €r 


a 


Bere. 5 
‘ be 54 -. - “a 
Pate 1 P (oa.2| 


‘ 


‘ a} 
c Fi vy 1] 
fae 8 OD | 
t') ix ‘ 


(0% relimke eagle 


be ; 1) vod Boge : 
2%) gnettcupe ialoneteae ts 


yee Wise éxsqmeay Le >. a 


37 


rates can be found and the rate of temperature changes can 
be reconstructed from the solution of oe differential 
equations, equations (2-5) are then defined if the increase 
in total temperature inside the combustion chamber due to 
fuel energy is known. 

Fig. 16 shows the schematic diagram of the combustion 


chamber. 
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Figure 16. Schematic diagram of the combustion chamber 


From the first law of thermodynamics: 


hic T + Woe T + h(t} = 210 
a5 Pa o5 ae Pe, 5 (t) ( )) 
CW yA sete W - W ya SPS Sern ham *)'L 
Pp fout o6 Brae fin fout o6 a a6é as) 06 

where: Cp = specific heat of air 
a 
Cc = specific heat of liquid fuel 
Cc = specific heat of burned fuel 
P£g 
h(t) = instantaneous fuel heating value as 


defined in equation (2-2) 
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Since T) and T. in equation (2-2) are very small 
quantities, it can be assumed that h(t) = Wh, where hy 
is the heating value of fuel per unit weight, since com- 
bustion chambers are usually designed to carry out complete 


combustion way up before turbine inlet. Thus, equation (2-10) 


in terms of fuel air ration W,/m reduces to: 


a5 
Wer 
p05 + Mas [c, Tor + he| 
TO6 = paar ce (2-11) 
ee Cc a7 (¢ 
a5 Prg Pa 


TO6 is a significant term since this is the highest temper- 
ature which the turbine inlet is subjected to. For control 
purposes, the measured TOG should account for: 

Ha transport delay , sede where s is the 
Laplacian operator and Ta is the time delay for the 
burning fuel to reach the turbine inlet, and 

2. time delay due to the inertia of thermocouples. 

It is worthwhile to mention here that among the 
quick acting means of temperature protection, fuel-air 
ratio limiters are widely used which are especially effective 
at starting and sharp load increase. However, fuel-air 
ratio limiters have low accuracy since changes in specific 


heat of fuel or changes in fuel gas ambient temperature and 


pressure do not make setting corrections. Thus there is 
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still a need for very quick acting thermocouples since the 

time constant for heating of the sharp trailing edge of 

turbine blades is in the order of one to three seconds. 
For the case of the LM2500 gas turbine, which has 

a single spool gas generator, equations (2-5a), (2-5b), 

(2-5g) and (2-5h) are simply deleted in the mathematical 


model. 
ITI-5 REDUCTION GEAR: 


As mentioned in Chapter I, reduction gear is the 
impedance matching component to assure that the power 
turbine rotor of the gas turbine and the pump rotor 
operate with minimum fuel consumption at a steady state 
on power levels most often used. Thus it can be represented 


in bond graph model as in Fig. 17. 
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Figure 17. Reduction gear bond graph model 
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where: Tp = power turbine torque 
@mp = angular speed of Sos turbine 
Tp = pump torque 
Wy = angular speed of pump 


m = modulo of the transformer TF or the 
reduction gear ratio 


I = combined polar moment of inertia of 
pump and power turbine referred to 
the pump shaft 


Re = shaft bearing frictional resistance 

Reg = reduction gear frictional resistance 
associated with reduction gear 
efficiency 


Assuming the reduction gear efficiency and shaft 


frictional resistance to be of the form: 


_W 
Cy 
= + & = 
oq = 0-75 + 0.22(1.0 - e ) (2-12a) 
= oF = as 
a Co C, WwW (2 12-b) 


and the constants Cy, C_ and C3 are to be determined from 


2 


data, then from Newton's Law: 


Ietges se imGeMad oxee 1% (2-13a) 
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For the case where the power turbine is developing 
negative torque at idle, fuel flow rate during transients 
equation (2-13a) should be: 


m Tm 
aden ae o (2-13b) 





II-6 PUMP: 


Pumps are energy converters transforming mechanical 
power to fluid power. In order to understand the basic 
mechanism of pumps, consider the velocity diagram repre- 
sented in Fig. 18 for axial and/or centrifugal pumps 
without guide vanes. The solid lines represent the design 
point velocity diagram while the dotted lines represent the 


off design velocity diagram. 


U = rotor blade velocity 


Wy = relative velocity of entrant fluid with 
respect to rotor 


c, = absolute velocity of entrant fluid to rotor 
and also exit fluid from stator 
Co = absolute velocity of exit fluid from rotor 


or entrant to stator 


W.. = relative velocity of exit fluid from rotor 
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Figure 18. Pump velocity diagram 


For a fixed geometry pump, referring to Fig. 18: 


aC, = W, (cos B - sin A tan7Y ) = O(B) (2-14a) 
Jak Ol fos Wi (cos B'- sin g'tan % = 2 (B') (2-14b) 
where: 
€ 
W, =u + ¢,° sin £ = ——— 
Vu2 + ei 
Parennee jc ; ~ e'5 
WwW i S sas I | sin g'= 
Uo es 


From the above equation it can be seen that the 
change in the tangential velocity component ( Ac.) of the 
fluid at any operating point is a strong function of 


dynamic similarity "flow coefficient" angle # and a weak 
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function of angle 8 since the deviation angle of the fluid 
at blade trailing edges is approximately constant. Since 
sin # and cosB are periodic functions, then pump and/or 
turbine characteristics must necessarily also be periodic 
functions of dynamic similarity "flow coefficient" angle # . 
From Euler's equation, torque (@) and pressure (P) 


are: 


T = paQrac, where r is the pump radius 


P= pr Aciw where W is angular speed 


Substituting equation (2-14) results in: 


i} 


© =[pr(p)] - 2=A(B) - 2 (2-15a) 


P=[er(p)] -P=-a,( 6) -o (2-15b) 


Thus, unlike positive displacement pumps and turbines, which 
are transformerly coupled, centrifugal and axial pumps and 
turbines are gyratorly coupled by Rg (AB Dire 

If appropriate power dissipation terms, Pa: are 


added to equations (2-15) such that: 


Ea ge Pozo (2-16) 


where: subscript pe= actual pump torque and pressure 


Pg power dissipation 
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«4 


then for a physical pump we have: 


Vas R, "Ot U5 = Bg NOt Rg w (2-17a) 
= = = = : = F: = 
Po “ Pa Rg w Ro Q (2-17b) 


dissipation torque 


Il 


where: aq 


= dissipation pressure 


tg 
Qu 
| 


The bond graph from representation cf equation (17) 
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Figure 19. Axial or centrifugal pump bond graph model 


where: Ll= polar moment inertia of rotors 


= fluid inertance in pipings 


H 
Fh 
| 


bearing frictional resistance 


a 


Rr pipe frictional resistance 


i} 


Note that the bearing frictional torque Rn is not the 


same as the pump torsional resistance R,,.- By referring to 
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4 


the typical centrifugal pump performance curves of Fig. 20 


gle 





Figure 20. Typical centrifugal pump characteristics 


and taking ncte of the points (x) where the power flow TW 
or PQ are zero across the gyrator (GY), it is evident 

: : d 2 
that the dissipation torque e g = T 3(Q, ) and the 
dissipation pressure a= (ona It can also be shown 
that Ry, and . are not independent and can be related by 
some function Rq(Q,W). There the bond graph of Fig. 19 


reduces to Fig. 21. 
r i 


eae, i 


% |e 


Figure 21. Axial or centrifugal pump model 
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where: R, = R(Q,00 ) = two port pump resistance 


bearing frictional resistance 


2 


w 
Il 


f frictional resistance on pipings 


From Newton's law and the above diagram: 


zit age = oe - Te. - (R,2 + RW) (2-18) 
go _ a: 3 = é 
I, ae (Rw Ro) Py Es (2-19) 


where: ts = pump input torque 


Cs = shaft bearing frictional torque 
Pee piping resistance friction pressure 
Po = pump pressure load 


For S.E.S. waterjet propulsion, the axial pump 
characteristic used in the system study (Fig. 22) was taken 
from ref. (6), where the head coefficient W = gip/U- is 
expressed as a function of flow coefficient f = c,/U. Here 


the pump pressure fe of equation (2-17b) is: 





2 Puzv 
H ern OG we) p=a sea RO! = RO (2-20) 
P g Pp g 3g Q 
where: p = density of fluid (pee) 


gravity constant (32.2 fey see") 


Q 
ll 











mo ietutua precre’ i mete 


teeipslh eveds sft tae wel a' ae 


(ef-s) “# = tga wa 


i ' 


mip2o3 deiet Spiteseg 


srotssia? pigseedt Jtade = 


G 


axlaeorves oboe s sx guandeiags praigiq 


taol stun hia gaa = 


cx Enizns e828 ,noretlmoig Jef toonw 


yas, SRM cf se Vode) Bees eet ext 


ry ang + ri th f me Seas mG: 
Ny - r2e200 © i. 3 ! 
* > 
f fats * 
7 P 
J 
a) 
pb 5 
Li 
we Se i 
= yi 


47 



















































































Ses Gxedd oad cased beets vase teeed cased ceeseeaees 

rit SEGEREEEEE 

\S8 SHES0 RSGEs eeeseeeees 

Seeertteet ceenteetns eeaeteeees 

pf 

Pence ences 

Sseesfeeeseressteent 

z Sa ee 

A sesey casas cnesi Gasaiatand Gass Gani Eon EEE 

oa fe Eetedirted Gerad eeed Cesetvenst 

hey ee eeeae neeene: age SUSUGSRGGGREGuGESRas 
oS 

aa Saasa feeetsezet 

fa See eeestael 


ad 











9 BSS88 seens eee 
ac Gooeceoeed cea 





ai 
aeaE 


















eogoem 
Soccdesoos 






























HEAD cOEFFICIENT YY 























Sines Sete snetaar enfant 
is ee ee SS eee soe ferereet eects eae 
Fes ieee afore or Gere teat toet 
EES EEN ie I 
Bebe fsen oleate 






fuiessssriiill 
led 


igs enes Spree cena 


FLOW COEFFICIENT 








FIGURE 22 PERFORMANCE GHARACTERIG TIC OF AKIAL 
PUMP 





a vee ™ wm soe vr 


48 


The design ¢ and Y of the above pump are 0.108 and 
0.51, respectively. To obtain ts of equation (2-17a), a 
dissipation function coefficient equivalent to equation (2-16) 


is assumed to be of the form: 


2 
= + - = 
eg = Ps. Cc, (¢- $,) (2-21) 
min 
where: Pg = dimensionless power dissipation 
[Ss ame ee minimum P . (at design point) 
min 


?3 = flow coefficient at design point 


C, = some constant to be determined 
relating to the efficiency charac- 
teristic of the pump 


Combining equations (2-21), (2-16) and (2-20), then 


Tp and P, of equations (2-17) in terms of pump dimension 


p 


reduces to: 


eS SAPO P es) an [257] * Dp (2-22) 
2 
a S ey leo y >| (2-22b) 
Be og 
where: ioe = pump torque (ft-lb) 

Y = dimensionless head coefficient 

g = dimensionless flow coefficient 

PGS dimensionless power dissipation of 


equation (2-21) 


P = water density (64 lb/ft?) 
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4 


A, = pump annulus area (£t2) 
= pump RPM 
D = pump blade tip to tip diameter (ft) 


Flushed and semiflushed inlets shown in Fig. 23 
and Fig. 24, respectively are considered to be more 
practicable to use for S.E.S. than ram inlets. 

The geometric parameters affecting the design 
pressure recovery and drag characteristics of flushed and 
semiflushed inlets are: entrance width to depth ratio, 
ramp wall divergence and ramp angle. For fixed geometry 
inlets, performance at various operating points is best 
described by efficiency of pressure recovery Ae asa 
function of velocity ratio V/V, where We is the velocity 
of the craft and V.. is the average axial velocity of the 
fluid across the inlet throat area A;. Some typical pressure 
recovery performances of flushed and semiflushed inlets 
(ref .(32)) are shown in Fig. 25 and Fig. 26. 


Similarly, the drag coefficient 
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RAMP ANGLE 


FIGURE 23 FLUSHED INLET 






STATION 2 


STATION 1 


FIGURE 24 SEMIFLUSHED INLET 
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FIGURE 25 FLUSHED INLET RAM RECOVERY RATIO FOR PARALLEL WALLS 
f - MEASURED AFTER INLET SECTION 
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FIGURE 2¢ SEMIFLUSHED INLET RAM RECOVERY RATIO 
I- MEASURED AFTER INLET SECTION 
2- MEASURED AFTER DIFFUSER SECTION 
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FIGURE 24 FLUSHED INLET DRAG 
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FIGURE 23 PRAG OF SEMIFLUSHEO IHLET 
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can also be expressed as a function of the ratio V/V, 
to describe the drag (D) Rip ee ee the inlet 
at various operating conditions. Some typical drag 
characteristics of flushed and semiflushed inlets (ref. (32)) 
are shown in Fig. 27 and Fig. 28. 

From the above pressure recovery characteristics 


the bond graph representation of inlet is: 


Figure 29. Flushed or semiflushed inlet bond graph model 


fluid inertance in the inlet section 


H 
i] 


where: 


£ 
Poi = total pressure of the fluid way 
ahead of the craft 
Re = equivalent resistance of the inlet 
plus the boundary layer losses at 
the sidewall of the S.E.S. in front 
of the inlet 
P42 = total pressure at the exit of the 


inlet or at the entrance of the 
diffuser section so that at the 


: = (le 
steady state P52 rel 


efnt aft 


: +5 
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II-8 DIFFUSER: 


The purpose of the diffuser is to increase the 
static pressure at pump inlet by converting the kinetic 
energy to static head. For a given expansion angle and 
fixed inlet condition, the total pressure loss of a 
diffuser is proportional to the theoretical value of the 


total pressure loss for a sudden expansion of the same 


Aout ) 


exit area to inlet area ratio (Ap = ae 


Thi - 
cP his pro 


portionality constant is called diffuser loss factor (Kg) 


and is defined by: 


poe Sines Feu) Sh (2-23) 
f (AR - 1) 
AR 
where: P| = total pressure at entrance to diffuser 
in 
= total pressure at exit to diffuser 
Pout 
q = dynamic pressure at entrance to 
diffuser 
De 
2p aoe 
qi 2g 
AR = area ratio 


Figure 30 shows the diffuser loss coefficient Ke as 
a function of area ratio (Ap) and diffusion angle (29) for 
Various geometrical shapes of diffusers. The bond graph 


representation, then, of the diffuser is: 
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RECTANGULAR DIFFUSER LOSS FACTOR, 
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Figure 31. Diffuser bond graph model 


so that at the steady state: 


Pa 1 : 
R 
Bae OS sate aS PY, 729 (2-24) 
II-9 NOZZLE: 


The purpose of the nozzle is to convert the excess 
static head of the total head by the pump into kinetic 
energy head so that thrust can be developed. A well-designed 
nozzle, unlike diffusers, has negligible loss so that 
efficiency of 96-99% can be easily attained. The bond 


graph representation of nozzle is: 
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Figure 32. Nozzle bond graph model 


so that at the steady state: 


ev; 
es a) 
2g = Mn (Po. - ma (2-25) 


where P., total pressure at the pump exit 


ain 


Pa = atmospheric pressure 


Un 


nozzle efficiency 


II-10 S.E.S., MASS, RESISTANCE AND WATERJET PROPULSION 

THRUST: 

The S.E.S. bare hull resistance data used in this 
study are tabulated in Appendix I. Note that if the 
boundary layer losses due to S.E.S. sidewall frictional 
resistance are included in the conputation for the total 
pressure recovered at the exit of the flushed or semiflushed 


inlet, then by energy conservation, an equivalent resistance 
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must be subtracted from the bare hull resistance data above. 
Otherwise, the propulsion system will eerecnatised twice. 
The assumed added mass for the forward motion of the S.E.S. 
was 10% of the total displacement. 

From the conservation of momentum, the instantaneous 


thrust developed for a fixed reference frame is: 


reece i [-& vav} I luietan,* a5] (2-26) 


or in terms of component thrust in x and y direction: 





pl, dQ Pe 

TY = Tg AL tge Ss cost - vo (2-29a) 
Pp L dQ P2 

ce = ar 3 ae + or (Vv. sin & ) (2-29b) 


where: L,. and Ly are the distances between inlet 
entrance and nozzle exit in x and y directions, respectively. 


At the steady state, the "power input" to develop 


thrust is: 





2 2 
PVo |V. 
2 TATE ere 4, -1]/f 0 (2-28a) 





i ee (2-28b) 
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where: 
APQ - TaV > 0 
a= 


Since 


then it is evident that the power conversion of equations 
Ovo lS 
(2-8) is gyratorly coupled with modulo m, = =a wes 1 


and the bond graph representation is: 


Is.e.s. 


ate ,—4 Rese ey i 1 
(Me 


Rg R, + Ry 


Figure 33. Bond graph model of thrust develop in S.E.S. 
waterjet system 


The equivalent dissipator resistance 
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where 
Vv V. 
Po ee 5 2 
AP off g WG - 1 ve ar a Ve 
P 2|V 2 
Pp, =+iy, E eal 
an 2 Gg 
oO 
is: 
ale 1 W5l_ V5 _ v; 1 
Ry = \vo| VY, eee - To - i Vi @} (2-29) 


The overall waterjet system bond graph model which 


combines all the different components is shown in Fig. 34. 
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CHAPTER III 
SIMULATION AND CONTROL OF S.E.S. WATERJET 


PROPULSION PLANT 
III-1 GENERAL: 


For the purpose of controlling the power turbine 
speed, the implementation of the detailed model of the gas 
turbine was found to be unnecessary. This.was primarily 
due to the higher frequency response of the combustion 
process and pressures in the storage volume of the gas 
turbine in equations (2-5) compared to the frequency response 
of power turbine speed. However, for a very low value of 
polar moment of inertia of the power turbine drive train as 
referred to the pump shaft, the poles and zeroes of 
equations (2-4) and (2-5) will be significant in the design 
of a stable high frequency controller. Likewise, if the 
primary concern of designing a controller is the avoidance 
of surge in the compressor and of overtemperature in the 
turbine inlet of the gas generator then the detailed model 
will be necessary. 

In the case study it was assumed that the LM2500 


gas turbine used in the simulation was adequately protected 
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from overtemperature and surge by its own control package. 
This gas generator protective control setup was simulated 
by putting a scheduled limiter on the rate on which fuel 
flow demand signal can be increased. This rate can be set 
to have lower value than the built in limiter in LM2500. 
The purpose of this limiter is to allow time for the gas 
generator rotor to build up speed and in turn air mass 
flow rate in order to neutralize the increase in pressure 
for surge avoidance and the increase in temperature due 

to sudden fuel flow rate increase. 

As for the rest of the gas turbine dynamics, it was 
noted that the instantaneous torque developed at the power 
turbine is a function of the instantaneous air mass flow 
rate and power turbine speed. Air mass flow rate in turn 
is dependent on the instantaneous temperature and pressure 
at the storage volume of the duct between power turbine and 
the last stage of the gas generator. Thus the transfer 
function of the power turbine torque from equations (2-4), 


(2-5), (2-7) and (2-9) can be approximated by: 


ss ve Tagstl Tyst1 Tae Tys2 + (TegtT,) stl * 


where: if a = instantaneous power turbine torque 


Rr 
n 
n 

! 


= steady state power turbine torque which 
is a function of fuel flow rate (We) 


os ; 
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a 


and power turbine speed (N,) 


a 


= time constant for the gas generator 


XR 
< 
I 


time constant for the volume between 
power turbine and gas generator 


Laplacian operator 


n 
i} 


Since ies and Ca care va Gre sete SF noe second, 
which is very much lower than the time constant of the 
power turbine drive train speed, the first term of equation 
(3-1) can be neglected so that: 


Tr al 
Tes Cag tty sri (3-2) 





The plot of ce as function of fuel flow rate We 


and the power turbine RPM is shown in Fig. 8. 
IITI-2 DYNAMIC EQUATIONS OF THE WATERJET PROPULSION PLANT: 


With the above approximation, the water propulsion 
system was reduced to the following non-linear ordinary 


differential equations: 


a. The S.E.S. Velocity Equation 


=a = [2,,(vy.9) See R, (V,,0,A4)] (3-3a) 
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where: V_ = S.E.S. velocity (ft/sec) 
W = S.E.S. mass plus 10% added mass (lbm) 
G.o= gravity (32.2 ft/sec’) 
T = waterjet propulsion thrust (lbf) 
R. = S.E.S. bare hull resistance (l1bf) 
R, = flushed inlet appendage drag (lbf) 


22 
A. = flushed inlet throat area (ft ) 


The instantaneous thrust Cr) in the longitudinal 


direction of the craft from equation (2-27a) is: 


pL dQ dv PQ 
y= —~= | — _ a +e nim. weOsalo= VV.) (3-3a-1) 
= g dt dt g 5 ° 


where: $ = water density (1bm/£t?) 


L,, = longitudinal distance between the flushed 
inlet entrance and nozzle area (ft) 


A, = effective cross sectional area in 
longitudinal direction (ft) 


V. = waterjet velocity at the nozzle exit 
(ft/sec) 


o& = depression angle of the nozzle 


The S.E.S. bare hull resistance (R,) should be a 
function of S.E.S. velocity (V, and S.E.S. operating draft 
Or air cushion pressure. In the simulation, however, it 
was assumed that the S.E.S. was at designed cushion pressure 


above hump speed. Below hump speed a resistance (R,) 
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i 2 
proportional to Vo was assumed, so that: 


m= 167.72 V2 0<V <50.7  (3-3a-2) 


2 
R.= 26.726 V - 6935.34 V_ + 712204.0 V 250.7 
s 12) e} Om 
The flushed inlet drag R from equation (2-23) is: 


eg 
R= hg 2 


0° ¥, (3-3a-3) 


For the purpose of controller design the flush inlet 
drag coefficient Cg was assumed to have a sharper point of 
minimum value at V3/Vo=0-8.- The minimum value of Cg was 
also assumed to be equal to 0.15 consistent with the published 


experimental data, so that: 


Cg = 1.0 V;/V,< 1.0 
Cqa = 1.0 - (0.8/0.85) (V4/V,) 1.0£V;/VK0.8 
Ca = 0.2 Vj/V, - 0.01 V;/V,7 0-8 


b. The Pump Speed Equation 


aN . 60 A i : 
dt amr, [ To(7 g-RoR) oe (N) v0.0) | (3-3b) 
where: ne pump RPM 


H 
iT) 


polar moment of inertia of power 
turbine, reduction gear and pump 
rotor referred to pump shaft (ft- eeees ) 
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@. = effective engine torque at pump shaft 
s 
(£t-lb£) 
Up = shaft bearing friction torque (ft-lbf) 


Tp = pump torque (ft-lbf) 


The effective engine torque at pump shaft from 


equation (2-13) is: 


Tv, = (RGR) ce ane T,20 (3-3b-1) 
e te “ 
s = (RGR) KO 


where: RGR = reduction gear ratio = 5.55 


MRG = reduction gear efficiency 

a -N, /895. 

‘MrRG = 0.75 + 0.22 (l-e ‘E ), where 
Ne = power turbine speed 


The shaft bearing torque (T,) is assumed to be: 


- = 1000 + 15N (£t-1bf) (3-3b-2) 


The pump torque (7p) from equation (2-22a) is: 
2 
PA, | N 3 
rv SAB pe Wee 
Coe e BA) = | pb (3-3b-3) 


where 


< 
n 
g = ID 
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A curve is fitted to the pump characteristics of 
Fig. 22 for the head coefficient (¥) as a function of 


flow coefficient (fp) so that: 


= dans Get - 1.916) + 0.85 PK0.04 
= 220 6° = 25.59°+ 1.308 0.0 <6<£0.065 
Y= 68.2176" + 10.173 + 0.2069 0.065¢$<0.130 


= -318.45¢7 4 72713 @ - 3.6949 0.1306$£0.152 


The dimensionless power dissipation Pe OF equation 
(3-3b-3) was normalized to a pump having an efficiency of 


90% at design point so that equation (2-21) yields: 


P, = 0.00623 + 12.4 (dt =sc tee) 


c. Waterjet Flow Rate Equation per Pump 


ets (0) 
See hy CNG +H V_ YA.) =H + H,. (—2 3- 
dat a a Q) eg “a 3) = (2) = Ab (3-3c) 
where: Q = water flow rate fee eee) 
3 
P = water density (lbm/ft ) 
‘ ; : ; 2 5 
I. = fluid inertance in pipe (lb-sec /ft ) 


Hgy = pump net positive suction head (ft) 
H, = waterjet kinetic head at the nozzle (ft) 


Hy, = head generated due to S.E.S. 
acceleration (ft) 
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«Aa 


The pump head Ale from equation (2-22b) is: 


Hy = . [eon p| : (3-3c-1) 


The net positive suction head Hoy from equation (2-25) 


and from flushed inlet efficiency characteristics is: 





2 
1 ¥, 
where: he flushed inlet efficiency 
Mp = diffuser efficiency 
Aur, = pump height above waterline 
-0.10/A;V, 
M4 = 0-85(1-e Oo) 
= Ap-l 2 
Qo = 1K) 


The waterjet kinetic head at the nozzle (H) xe: 





2 
Ae 
where: Mn = nozzle efficiency ~ 1.0 
V5 = Q/C_Anoz = waterjet velocity (ft/sec) 
ee nozzle coefficient of area con- 


traction ~ 1.0 


Ayoz = cross selena area of nozzle 
throat (£t”) 
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The head generated through S.E.S. acceleration (Hy) 


is: 
H zs fae ee (3-3c-4) 
— As -3c- 
Vv g | at A, dt 
where: L = distance between entrance to flushed 


inlet and nozzle exit in the we 
direction (ft) 


d. Power Turbine Torgue Equation 


— #_ 9.5 | Es MerNy) - ies) (3-3d) 


= 
(0) 
K 
: 
R 
Il 


instantaneous power turbine torque 
(ft-lbf) 


st 
I 


= steady state power turbine torque 
as found in the engine performance 
map 


N = power turbine RPM 


Mie factor 0.5 is the iS Sees) term in equation (3-2). 


e. Fuel Flow Rate Equation 


aw 
ae =r2e0 


Fh 


(Wea - W,) (3-3e) 


instantaneous fuel flow rate (lbm/hr) 


Sy 
0) 
K 
>) 
= 
i] 


= fuel flow rate demand signal (lbm/hr) 


Fis 
Q 
| 
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Equation (3-3e) assumes that LM2500 control package 
has a built-in servomechanism for the fuel pump with a 
time constant of 0.5 second and a constant gain of 1.0. The 
rate limit and hard limit were incorporated in the fuel flow 
rate demand signal so that 
W..<W 


cae eae 


<W <w 
foun 7 fa ida. 


The minimum fuel flow rate demand signal Wea was 


set to idling fuel flow rate of 800 lbm/hr and the maximum 


fuel flow rate demand signal Weg was set to 10,000 lbm/hr. 
max 
The maximum rate of change of the fuel demand signal Wr, 
max 
was set to six seconds from W. to W : 
fdnin fdmax 


£f. The Flushed Inlet Area Servomechanism Equation 


dA; 
“ar = 9-5 (Aig - Ay) (3-35) 
where: A; = instantaneous flushed inlet throat 
area (Ge) 
aae = flushed inlet throat area demand (ee) 
i 


Equation (3-3f) assumes a hydraulic servoactuator for 
the flushed inlet throat area with a time constant of two 


seconds and a constant gain of 1.0. In the simulation, the 
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hydraulic actuator was set to have a hard limit of 7.0 £t? 


> . 
maximum, 1.5 ft minimum, and a rate limit of 0.5 £t> /sde 


so that 

A <-0).5 
LS: A5 = 7.0 

III-3 PUMP CAVITATION: 


To simulate pump cavitation during transient 
condition at various pump operating points, a dimensionless 


suction head criteron was assumed to be: 


where: H = critical net positive pump suction 
SVC 
head (ft) 


Gq 
ll 


pump blade tip velocity (ft/sec) 


A curve is fitted into the pump cavitation data 


ig. 22) so that Y suc is defined as a function p - The 


pump is assumed to cavitate if the actual dimensionless 


suction head VES Bo Mey “fer below Wow: The degree 
u 2 


of degradation in pump performance was assumed to be of 


the form: 


SVC 


Ysvo-Ysva] 
ae = 2 c) (Sy) on Wee (3-5) 
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73 
where: Yeo = dimensionless head with cavitation 


y= dimensionless head without 
cavitation 


C, = some constant to fit the pump 
cavitation data 


Substituting equation (3-5) to equation (3-3c-1), 


the pump head generated with cavitation is: 


¥ Vsve eae N : 
i, = ae ey Se aa >| Ysva <¥sve 


(3-6) 
The dimensionless power dissipated when there is 
cavitation Pac is also assumed to be of the form: 
RP, =P, {l.0+C ——s baa ae (3-7) 
ed 2 eee SVA <TSvc 
where: Py = dimensionless power dissipated without 
cavitation as defined in equation (2-21) 
Cc. = some constant to be determined from 


pump data 


Substituting > and P in equation (3-3b-3), pump 


dc 


torque with cavitation is: 
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III-4 LINEARIZED DYNAMIC EQUATION OF THE WATERJET PROPULSION: 


The linear expansion of equation (3-3a) about a 


steady state operating point yields: 


= [Se T OR, 














ORG OR. OR, 
= av WeVed = ooo (Q-9.) = OA (a,-A;0) | 
The subscript 5 represents the steady state operating 
point values of the state variables and the partial deri- 
vatives are evaluated at the corresponding steady state 
operating point values. Similar expressions can be derived 


for the rest of the equations of (3-3). At steady state: 


Defining: 
e e e e 
Bem ac ‘ox ae Von cern os = os 
Q-2, = om and Q-Q, a or =O 


etc. yields the linear form of equations (3-3) 
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The block diagram of equation (3-8) is shown as 
solid lines in Fig. 35. The broken lines represent the 
eoncroller. The cuts in the diagram at points x, y, 


and z will be discussed in the succeeding sections. 
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III-5 FLUSHED INLET THROAT AREA CONTROL: 


For the performance of semiflushed or flushed inlet 
to be optimum 1) pressure recovery must be maximum and 
2) drag must be minimum. Fortunately, these two criteria 
were not contradictory, as can be seen in the performance 
Burves of Fig. 25, 26, 27, 28 .of ‘Chapter IT. At ithe inlet 
to ship velocity ratio of V;/V, = 0.8, the above two con- 
ditions are satisfied. A control algorithm then can be 
set up in such a way that V7" (or equivalently 
UAV, = 0.8) at steady state and at least close to 0.8 
during transients. A preliminary proportional and integral 


control for the flushed inlet actuator was set up in the 








form of: 
mem, Ce ( - 0.8) (3-9a) 
j = Be s-0:.8) (3-9b) 
idy In ‘AiV, 
so that: 
— ‘ + Ss 
ia Aidp Pe ae (3-9c) 


Substituting (3-9) to (3-3f) yields: 





e 
A; = -0.5A; + 0.5Aiq, + pees. (ASV, - 0.8) (3-10) 
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Taking the linear form of equations (3-9b) and 


(3-10) and including them in equations (3-8) yields: 











v, Bij 20 Aya, Of 0) ay, 0 [A y, 0 (3-18) 
Ny GE Raa: Ao 0 Ol EN 0 
Q, arenes 10! Bd Dh OO Q, 0 
rent fo) ig ce ea Ay. Os) Ae Cie | #10 [Fea | 
Wee O Bow OiVyD As. 0 ' © Wey hey 
Ais Ae, 9 Ags 0 0 Ace Aer] [Ais 
sa,, a 0 Be 0 8 Ae 0 Aids. 0 
where: 
33 aa,v| Q 
Ap = 0.5 Kon 3a), = -0-5 ca Poel 
Ore 
a[asva| sl 
pone OVSiKE) PeSegoee* hinds "palais 
; oe 
A. o= Aggt0.5Kpa—a5—> = -0.5 eal 2 |+ Aes 
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This PI control is shown in Fig. 35, in which: 


mel 
= FEES 


To obtain the steady state values of V ,N,Q, 
oo ° fo) 
eS ; hes and Ba an open loop response to different 
levels of Weg using conservative values of K and K 
IA PA 
were obtained. Some typical responses are shown in Fig. 36. 
To obtain an optimum Kia and Hat an open loop 


transfer function Ai4/Vo by cuttings the.loop at point. x 


* 


in Fig. 35 was obtained using MBODE program. Such open 


loop transfer functions can be represented in matrix form: 
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V, cs, 8 OO. Ge Vv 1 
N, Opes Ao 5 AG, 0. 00°10 NU 0 
Qo, A3, A, Ag3 9 29 0 Q. 0 
te} =| 0 Aaa 9 Aggy Ags 0 OO] | Tz, | +] 0 Ga 
a Gye 00h Ree Oe Ol are 0 
Ay Ag, 9 Fez 9 O AgeAga] | A, 0 
Aiat. Btn a) DG ee 0 
ly] = [o Gm arGe 10, 60.) 1 0| ve 

Ny 

Q, 

(ae 

We: 

ee 

idy 


where Aig of equation (3-11) is set to zero. 

The zero location can be adjusted by the ratio of 
Kya/Kp a depending on the location of poles and zeroes of 
the hydraulic servoactuator for the flushed inlet at various 
operating points. The hard limit and rate limit of the 
actuator on the other hand, sets a constraint on the 


Maximum value of Kra that can be used, otherwise excessive 
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overshoot occurs. For the case study presented, the hard 
limit and rate limit on the actuator mee teveecome by 
setting Kya = 0.5 and K, = 5.0 and putting a hard limit 
on the integrator portion of the controller equal to the 
hard limit of the actuator, so that the maximum and 
minimum values of “at were 7.0 and 1.5, respectively. 
Thus for the control of the flushed inlet geometry 
some means to measure Q, V, A; must be provided in order 
to keep track of V;/V,°- The time delays associated with 
these measurements must also be included for the proper 


location of the zero of the PI controller. 
III-6 OPEN LOOP SIMULATION TO FUEL DEMAND SIGNAL: 


The fuel demand scheduled limiter used in open loop 
simulation was set to have a maximum rate of change from 


idling Weg to maximum W in six seconds and a hard limit 





fd 
of idling Weg = 800 lbm/hr and maximum Weg = 10,000 lbm/hr. 
: : : Pox dQ Os 
The Ey term in the flow equation (3-3c) and the g |at-Ax dt 


term in the thrust equation (3-3a-1) were neglected since 
they were negligibly small compared to the other terms. 
Since S.E.S. under study uses six engines, the thrust 
equation (3-3a-1) and the inlet appendage drag equation 


(3-3a-3) were multiplied by six for the simulation of the 
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straight line motion of the craft. Computer runs were 
made using a stepped fuel demand signal (Weg) from each 
engine idling fuel flow rate of 800 lbm/hr to various 
levels of Wy, up to maximum We of 10,000 lbm/hr per 
engine. The stepped Weg Signal was filtered by the maxi- 
mum rate change of the Scheduled ee limiter with the 
setting discussed above. For a realistic open loop simu- 
lation, the flushed inlet area control loop discussed in 
the preceeding section was closed. Pump cavitation was 
observed at the stepped Ls of 5,000 lbm/hr/engine and 
above. Some responses of pump RPM, pump flow rate (£t2/sec), 
S.2.S. velocity (ft/sec) and flushed inlet area (£2) to 
the stepped Wr, made are plotted an Fig. 36% 

Fig. 36 responses, however, were obtained using a 
higher value of fluid inertance (I,) in the pipings than 
the computed value for a faster YSYS solution (i.e., using 
higher value of time step). Comparison of the response with 
the actual value of I (i.e., using very small value of 
time step) yields a similar result. The higher value of 
I, used served as a filter to the oscillatory solution of 

YSYS to pump flow rate during transient and especially in 
the cavitation regime when the time step is greater than 


0.1 second. 
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III-7 PUMP RPM CONTROL: 


Since pump shaft RPM appears to be the easiest to 
measure accurately and appears to be the most effective 
indicator of power turbine stress overload, it was decided 
to base the fuel demand signal on the error (Ey) between 
the throttle RPM setting (Np) and the actual pump shaft 
RPM (N). 

By inspection of Fig. 35, the open loop transfer 
function Ny/Wegy yields a type O system. Such open loop 
transfer function, represented in matrix form of equation 
(3-11) with the output matrix [c] = [0100000], was 
verified using MBODE program and the result yields no pole 
at the origin in the s plane. Thus, to reduce the error 
(Eyx = Ng N,) to zero at steady state, an integral 
control is needed. A proportional control is also needed 
to bring N, close to Nr, faster for large error (Ey) - In 
order to obtain the gains Kone Kin for above PI control, 
MBODE program was used to locate the open loop poles and 


zeroes of N,y/Weg, for various levels of W i.e., 


fd' 
800 lbm/hr to 10,000 lbm/hr. 
In the MBODE solution, however, of equation (3-11) 
i A 
the A, term was so big compared to the Aly Ayo: Aaat Be 


and Age terms that a big numerical error resulted. It 
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was therefore necessary to reduce the order of the system 


to six by setting: 


This is equivalent to setting the fluid inertance to zero 


or shifting the pole due A to -—e9 in the s-plane. Such 


33 
a procedure will yield a solution very close to the true 
solution since the pole due to aoe being so far away to 


the left of the other system poles and zeroes, is effect- 


ively at -c&®. Solving for Q yields: 
* 5 


A A 
o, =v, -—-2N (3-13) 
A A = 
33 33S) 
Substituting equation (3-13) into Vv P , Ae and A. 
* Ex a id; 
of equation (3-11) results in: 
Vs a4 <5 0 6) ane 0) We 0 
N, Aa51 Ao a53 O O O Nv O 
Tt 0) a a33 a 0 ) 96 0 
Be 32 34 a ca 
Me ) 0 Gy edhe O}| We by 
an Se She i Ma Ma: (ig (ess oe 
A, 4 ao. 362 © b> a © A 0 
I I, 
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With the above PI control, since W = 2.0 (We 5-W.) 

K +E 
IN 

and Weg = (oy : eae), (N-N), then the closed loop 





transfer function N,/N, is represented by: 
* 


wi a Oe 0. ay, 00 Vv, 0 (3-15) 
Ny Q51 855 423 O O O 0 N, O 

Te, 0 a5 8ag e549 «0 C0 ou 0 

ie. =| 0) ay, a0% a, 0 0 a,, We ey I] 
As. Ac 452 0 O a55 S56 O ae 0 
ar, Se tae? A. 0. 0 P idye 0 

fd,, Glare oO 0) 0) 6 Wray Ee 


with the output matrix {c] =[o 10.0: 0.0 0] 


where: 


ao = -2.0 Kpy 
an ee) 
S72. TN 

7 = +2.0 Koy 
b = +K 
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An MBODE result for the value of matrix A of 
equation (3-14) at Weg = 5,000 lbm/hr with output matrix 
[c] = [0 1000 0] indicates that three of the zeroes 
cancel out three of the poles of the open loop transfer 


N (s)/w (s) and is plotted in Fig. 37. 
# fd, 
je 


-2.0 TOOTS Cone rom 





Figure 37. Open loop poles and zeroes of N,/Wea, 
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Thus in the PI controller of the form: 


Kin = Koy | Kr 


a aa S pn ‘ 


Kin/Koy can be adjusted so that the zero of the controller 


is in the vicinity of -0.234 and the gain K 


PN also adjusted 


so that closed loop pole is within the 45° line indicated 


in Fig. 38 for a fast response with less than 5% overshoot. 


: e 


aera * 
; FOR 5%. 
/  ovensnoot “\ 
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Figure 38. Root locus plot for K K.. =. 0.234 
Gg Pp tw PN 
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III-8 ANTIWINDUP LOGIC FOR PUMP RPM CONTROL: 


In order to compensate for 1) the excessive over- 
shoot and oscillation during transient caused by the hard 
limit and rate limit on Weg signal and 2) the shifting of 
the two rightmost uncancelled poles of equation (3-14) 
towards the origin at low power level, the Sor and Kix 
were adjusted accordingly in the form as shown in Fig. 39. 

The idea behind such a scheduled Kyy and Kpy isto 
set the gains of the PI control for small error signals to 


be equal to the optimum values as found in the linear 


analysis of section III-7, so that: 


> 
I 


£, (N,) 


ya 
It 


ola 


The £, (N,) and £, (Np) are the plots of the optimum ae 

and Kin as function of N. which will yield a fast stable 
response with less than 5% overshoot about an operating 
point 7 for small error disturbance. For large error 
disturbance, however, the value of K ah is too large to 
be compatible with the rate and hard limit of the scheduled 


fuel demand limiter. Large overshoot occurs because by the 


time the integral part of the fuel demand signal (Weg ) has 
I 
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Key = KomMay (1-0 ~ ABS (Ne-N)/100-0) 


ABS CNe-N) 


Kin 


Kru = Krwax (1.0 - 48S C(Ne-N) /100-0) 


ABS (Nea-N) 


FIGURE 39 PpuMe RPM PROPORTIONAL AND INTEGRAL 
SCHEDULED GAINS 
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reached the correct value of We for a steady state Np. 


the error signal is still lagging too far behind. This 


is due to the fact that the correcting signal Weg acting 


on the fuel pump servo is considerably less than the 
combined integral and proportional fuel demand signal 


Weg ) due to the filtering action of the scheduled 
Pp 


Weg limiter. Since for large error the proportional part 


(W. + 
fd; 


of the fuel demand signal (Weg ) by itself is large enough 
P 


to be filtered by the W 


fa limiter, K therefore can be 


IN 
set at a lower value and still get the same speed of 
response in pump RPM (N). Thus K, , can be set so that 
min 

by the time the error signal goes to zero, the integral 
part of fuel demand signal (Weg) is equal to the actual 
fuel flow rate (W.) for reference pump RPM (N,) - Over- 
shoot through this setup is eliminated. 

However, for Kp = constant = K , it was observed 

max 

that some secondary oscillation occurs during the accelera- 
tion transient (Fig. 33b). This is due to the hard limit 
on the Weg signal for large error. When the error signal 


was decreased to a value such that Wea was no longer hard 


limited, the Weg Signal decreased rapidly, which caused 


the actual fuel flow rate We to decrease slightly. To 















‘Sahoterie air to aibtos peisettas eta: fa oud a 
i in 7 , 

esq Innehaeoqetg ads roves wptst eH ants : 

deuce epret el iisest yi t seegith lavsle’ Samed 


Ce 
_ , 9 
a 
os 


a 7 


se nse Oiotored? 4 TF  rayint) )." oe ye 708 


*> bepoe smsc ef? fap fLi36 Sea oorlov 


es 


eM oa jax ot nae Ss ay! .(e) Me qmg 
; eta wd 


(axpoGni off .o1e2 oF ségh [Anyis weave off: 


fastos alis of Latpea ef (_eaw) fénpie Srtheb 5 


sav0 (-4) WEE qgietvg ooneuele: ick |.%). ede ag 


~ > 
Hoaserkagiea af quiae «lid mp 
: aes | A a, ' - Faas b | 
mt 
_O-S | » rise hI 71 $y Caste 
4 it fk aauge 
a oper! Leavin Ww 
: 53 
TT ; c of hemetrosm 
£ t+ bBeator>: xe w odd .boedl 


ip 


> Aa \ c.g? vi 5 Ww lou? innijon SOF) 


98 


overcome this effect, a negative SUS TS RS of Kay 
with respect to absolute value of error was set so as. to 
maintain always an increasing fuel flow rate (W,) during 
acceleration transient. ert was also set so that for 


large errors Weg will still be hard limited. 
P 


III-9 CONTROL ALGORITHMS TEST RESULT: 


The following hard maneuvers were simulated to test 
the control algorithms for pump RPM and for the flushed 


inlet throat area: 


(a) Acceleration from 109 to 635 RPM Throttle Setting 


This simulation was conducted by putting a step change 
in RPM throttle setting. The 109 RPM corresponds to idling 
power setting and the 635 RPM corresponds to 95% of the 
maximum intermittent power rating of the engine. The results 
are plotted in real time (seconds) in Fig. 40. The labels 
and scales are: 

V ="S.BUS. velocity (£t/see) 

N = pump RPM 

Q = water flow rate per pump (FE Jace) 

Number of pumps = 6 
WF = actual fuel flow rate per engine (lbm/hr)X100.0 


Number of engines = 6 
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AI = actual flushed inlet throat area per inlet 
(£t?)xX50.0 | 


Number of inlets = 6 


TRST 


total thrust of the six (6) waterjet system (lbf) 


TE = engine torque (ft-lbf) 


TP pump torque (ft-lbf) 


Notice that the gas turbine is dually protected from 
surge, overtemperature and transient mechanical overstresses 
by 1) the fuel demand rate limiter, and 2) by the PI con- 
troller gains scheduler of Fig. 39. Thus for large errors 
the fuel demand signal received by the engine will at most 
be a six second ramp signal. This is verified by the WF 
plots in Fig. 40. 

The response in pump RPM shows that the controller 
setup of Fig. 39 is fast and stable. In addition, it has 
no overshoot. There is also no secondary oscillation due 
to the rate and hard limits of the Weg scheduler during 
transient as this was taken care of by the negative sloping 
function of Kon with respect to the absolute value of error. 
Thus as the error decreases the fuel pump still receives an 
increasing Weg signal. 

The flushed inlet throat area oe increased in the 


initial transient to accommodate more flow at low S.E.S. 
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a 


speed to generate the acceleration force. It finally 
decreases as the S.E.S. builds up speed, maintaining NGS Y 
close to 0.8 value during transient and exactly 0.8 at steady 
state. Thus, setting a hard limit alone on the integral 
part of the PI controller for the flushed inlet is enough 
to overcome the rate and hard limit of the flushed inlet 
actuator. 

The pump torque (7) and engine torque ( c.) are 
plotted to check mechanical stress overloads during 


transients. 


(b) Deceleration from 635 to 150 RPM Throttle Setting 


This simulation was also conducted by putting a step 
change in RPM setting. In order to protect the gas turbine 
gas generator turbine inlet from sudden change in tempera- 
ture the fuel demand rate limiter for decreasing fuel flow 
rate was set arbitrarily to six second samp. The PI con- 
troller gains scheduler of Fig. 39 likewise contributes to 
this protection as an added precaution. The results are 
plotted in Fig. 41. The labels and scales are the same as 
those of the acceleration test. 

From the plots it can be seen that the response of 
pump RPM was slower than in acceleration and has slight 


overshoot. This was expected since the damping term in 
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pump RPM varies approximately to the square of pump RPM. 
Such an overshoot, however, is within acceptable limits 
and in addition it occurs at approximately 10% of the 
maximum intermittent power rating of the engine. 

Thus, it is evident from the pump RPM responses 
that the PI controller gain of Fig. 39 works optimally 
throughout the whole operating range of the waterjet system, 
that is, 1) fast response, 2) less than 5% overshoot, 

3) stable and 4) no steady state error. 

From the plot of We it can be seen that the actual 
fuel flow is rate and hard limited even without the actuator 
itself being rate and hard limited. 

The plot of the flushed inlet throat area (A, ) showed 
an initial decrease to its minimum hard limit for several 
seconds in an attempt to bring WIV equal to 0.8. This did 
not present a problem on the integral part of the PI con- 
troller for the flushed inlet since the integrator was also 
set to have a minimum hard limit equal to the minimum hard 
limit of the flushed inlet actuator. As the water flow rate 
Q and S.E.S. velocity Vag came to a steady state, the inlet 
throat area (A;) increased to its steady state value so that 
WitVs is exactly equal to 0.8. 

The pump torque (T) and engine torque (ze were 


plotted to check back flow of power during deceleration. 
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CONCLUSIONS 


From the result of the open loop simulation to a 
step fuel flow rate for the waterjet propulsion for S.E.S. 
several conclusions can be made. The first conclusion is 
that despite the pump being designed not to cavitate at 
hump speed with a thrust to resistance ratio of 1.2, cavi- 
tation can still occur during transients since designed 
pump head and designed pump flow rates for the hump speed 
can be reached way earlier before the S.E.S. goes to hump 
speed during acceleration. However, such cavitation only 
occurs for a few seconds. 

The second conclusion is that for a fixed designed 
area nozzle the pump is operating very close to the designed 
qd and WY at steady state for all power levels; likewise, 
the gas turbine is operating very close to the point of 
maximum efficiency for all power levels at steady state. 

The third conclusion is that in the criteria for 
designing the speed at which the inlet area can change, there 
are two variables involved, namely the flow rate Q and the 
S.E.S. velocity Vo° During transients it was noted that Q 


changes faster than V,. Thus, if local cavitation is to 
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be avoided in the inlet during transient, the rate of 
change of the term M/AV, should be considered. 

As already discussed in the MBODE analysis of poles 
and zeroes locations of the system the prime mover and the 
pump dynamics can be decoupled from the S.E.S. motion 
dynamics so that the design of a stable control for pump 
RPM can be simplified. A similar finding was found to be 
true also between the inlet area dynamics and the S.E.S. 
motion pease 

In the fuel demand scheduled rate limiter cebeacbeesied 
in the RPM ae of pump, cavitation can be avoided during 
transient by setting the rate limiter to some lower value 
for the normal operation of the S.E.S. For the hard accel- 
eration maneuvers the fuel demand scheduled rate limiter 
Gein be set according to what the gas turbine characteristic 
may allow for surge and overtemperature. 

The control algorithm of Fig. 39 works for a rate 
limited and hard limited actuator and is recommendable 
since analog or digital implementation is relatively simple. 
It also works throughout the whole nonlinear operating range 


of the waterjet system. 
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RECOMMENDATIONS 


The paper presented as it is, is like a textbook 
approach to the basic principles involved in the system. 
Several engineering approximations were made for a sim- 
‘plified digital simulation and control of the pump RPM and 
flushed inlet area. The following recommendations are 
then made for future work in this area of study. 

1. In the control of the flushed inlet, anes the 
throat area was considered in this work. This was in- 
tended to serve as a reference in the actual reshaping of 
the whole flushed inlet geometry consistent with maintaining 
maximum possible pressure recovery, avoidance of cavitation 
anid minimum drag. The mechanical implementation of this 
scheme deserves some further study. 

2. The fuel demand scheduled rate limiter still 
needs some refinement to avoid cavitation during acceleration 
for S.E.S. normal operation. Such a limiter should be 
consistent with the operating height of the pump above the 
waterline, and S.E.S. motion, prime mover and pump dynamics. 

3. In order for the S.E.S. to maintain a straight 


line motion there is a need for another set of controls to 
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regulate the thrusts of all the waterjet systems. Similarly, 
turning maneuvers with and without the use of thrust de- 
flectors deserve some further study. 

4. Backing maneuvers with thrust reversers were 
not included in the study due to the absence of data on 
flushed inlet performance in the backward motion. This 
problem should also be investigated. 

5. It will be recommendable to implement a hybrid 
simulation of the system so that the complete gas turbine 


model as presented in Chapter II can also be included. 
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APPENDIX I 


The following are the particulars of the S.E.S. 


under study: 


Length to beam ratio = 2.37 

Cushion pressure to length ratio = 2.5 
Cushion pressure = 412 lbt/£t- 

Displacement = 2109 tons 

Number of lift engines = 2 

Number of propulsion engines = 6 

Designed cruise speed = 80 knots 

Maximum intermittent engine rating = 22,500 HP 
Maximum continuous engine rating = 20,500 HP 
Range = 2,000 miles 

Hump speed = 30 knots 


Bare hull resistance (smooth sea, cushion borne): 


Velocity Resistance Velocity Resistance 
(Knots) (1lbf£) (Knots) (1bE£) 
28.45 422918 45 299578 
30 429438 50 282030 
35 367391 55 271360 
40 326606 60 265899 
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oa 


Velocity Resistance Velocity Resistance 
(Knots) (1bf) (Knots) (1lbf) 
65 261984 85 269772 
70 258695 90 278874 
75 259040 95 290098 
80 263038 100 303252 


Geoeds 


#entas 
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APPENDIX II 


The steady operating point values of S.E.S. velocity 
Si pump RPM (No). water flow rate per pump (Q.4). 
throat area per inlet (A; ) and reduction gear efficiency 
fo) 


Os ) for different values of fuel flow rate per engine 
° 


(We) are: 

oe Veo Ny Q, Aio m 
(1bm/hr) (£t/sec) (RPM) (et?/sec)  (£t7) RG 

800 SF 109 47,38 6.307 0.858 
2,000 20.1 241.8 101.6 6.318 0.92 
3,000 28.96 346.8 146.4 6.319 0.944 
4,000 34 a7 410.3 173.7 6.317 0.952 
5,000 43.22 465 .6 a i Basis 0.957 
6,000 46.62 514 218.0 6.317 0.961 
7,000 49 .67 553.9 235.6 6.317 0.962 
8,000 103.6 589.9 2522 6.317 0.963 
8,500 124.0 605.8 263.9 3.202 0.964 
9,000 124.0 620.9 280.4 22822 0.965 
9,500 1£35..2 635.6 292.0 2.002 0.966 


10,000 142 647 .9 296.6 2.628 0.966 
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a 


From the above operating point values, power turbine 


torque (T, ), pump flow coefficient (B.) and head coef- 
° 


ficient (yw ) can be determined. Thus the matrix A _ of the 
° 


linearized waterjet system equation (3-8) can be evaluated 


so that: 


-0.0264 


for Weg = 5,000 


-0.1088 


0.0 


-0.0293 


1225 


0.0 


-0.0106 


Swe 


-20.179 


0.00514 


-0.00757 


-3.024 


0.0 


0.02139 


-0 .04479 


-12 .337 


0.0 


0.0 0.0 -0.00055 


0.001302 0.0 0.0 


0.0 0.0 0.0 
-0.5 4.0 0.0 
0.0 -2.0 0.0 
0.0 0.0 -0.5 
0.0 0.0 -0.0094 


0.001449 0.0 0.0 


0.0 0.0 0.0 
-0.5 3.4 0.0 
0.0 -2.0 0.0 
0.0 0.0 -0.5 
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for W__ = 10,000 
-0.03283 0.0 0.02556 0.0 0.0 -0.1242 
0.0 -0.1754 0.0019 0.0014 0.0 0.0 
ie 7.704 -19.706 0.0 0.0 0.0 
0.0 ~32.506 0.0 -0.5 2.65 0.0 
0.0 0.0 0.0 0.0 -2.0 0.0 
0.0 0.0 0.0 0.0 0.0 -0.5 


The above constants were used in designing controller 


gains for small disturbance. 
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APPENDIX III 


Subroutine EQSIM is required by DYSYS for the solution 
of the dynamic equation of the system. The DYSYS program is 
available at M.I.T. mechanical and civil engineering computer 
facilities. 

Subroutine LM250 represents the LM2500 performance 
map of Fig. 8. The inputs to the subroutines are power 
turbine RPM and fuel flow rate and the output is power 


turbine torque. 
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